Quantum Dynamics of the Hubbard-Holstein Model in Equilibrium and 
Non-Equilibrium: Application to Pump-Probe Phenomena 
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The spectral response and physical features of the 2D Hubbard-Holstein model are calculated 
both in equilibrium at zero and low chemical dopings, and after an ultra short powerful light pulse, 
in undoped systems. At equilibrium and at strong charge-lattice couplings, the optical conductivity 
reveals a 3-peak structure in agreement with experimental observations. After an ultra short pulse 
and at nonzero electron-phonon interaction, phonon and spin subsystems oscillate with the phonon 
period T p h ~ 80 fs. The decay time of the phonon oscillations is about 150 — 200 fs, similar to the 
relaxation time of the charge system. We propose a criterion for observing these oscillations in high 
T c compounds: the time span of the pump light pulse r pumv has to be shorter than the phonon 
oscillation period T p h- 

PACS numbers: 71.38.-k, 71.10.Fd, 72.10.Di 



It is accepted that the physics of the high critical 
temperature (high T c ) cupratc superconductors is con- 
trolled by the interplay of different kinds of elementary 
excitations. The Mott-Hubbard nature of these systems 
was quickly recognized, and highlighted the relevance 
of the electron-electron interaction (EEI) [l[ . Angle re- 
solved photoemission spectroscopy (ARPES) Q and op- 
tical conductivity (OC) Q measurements revealed the 
importance of the electron-phonon interaction 

In spite of the general consensus that both EEI and 
EPI manifest themselves in the physical features of these 
compounds, it is still unclear which of these two in- 
teractions is responsible for particular properties. This 
difficulty arises because the EEI and EPI have similar 
spectroscopic energy scales. Ultra-fast pump-probe spec- 
troscopy, which probes the time dependence of the prop- 
erties of a system after it was displaced from equilibrium 
by a laser pump pulse, is a promising method for disen- 
tangling charge, magnetic, and lattice degrees of freedom. 
Indeed, since laser pulses couple to charge excitations, 
the characteristic response times of the different subsys- 
tems are not determined solely by their own characteris- 
tic energies but are strongly influenced by the strength 
of their entanglement with the perturbed charge density. 

Although pum ping in pump-probe experiments on high 
T c compounds [5l4ll| or any other material (see, e.g., or- 
ganic charge-transfer crystals 12-14j]) is always done by 
an ultra-short pulse of light, there is a significant varia- 
tion in the experimental setup. Conditions that can vary 
include fluence, photon density (total number of pho- 
tons per site absorbed during the pulse) , frequency of the 
pump photons 0J pU mp, pulse duration T pump , time resolu- 



after the time delay. Examples of such quantities are lat- 
tice structure and dynamics measured by ultra-fast elec- 
tron diffraction 0, Q and electron temperature derived 
from ARPES [(§. However, one of the most informative 
pump-probe experimental techniques is associated with 
energy- and time-resolved analysis of the electrodynamic 
response, such as reflectivity @, M, EMI , 

transmission 

HHl, or full OC In particular pump-probe experi- 
ments on high T c compounds [5l4llj revealed three basic 
timescales. The charge immediately reacts at times < 50 
fs. The phonons which are directly coupled to charges 
are involved in time scale of 100-200 fs, and the energy 
is gradually dispensed to the whole phonon subsystem 
(heat reservoirs) by anharmonic interactions during fur- 
ther 500-1000 fs. 

In this letter we study the OC of the two dimensional 
(2D) Hubbard-Holstein model (HHM) both in equilib- 
rium and after an ultra-short powerful light pulse. We 
choose this model because it reproduces the equilibrium 
OC of the doped compounds, where experimental stud- 
ies resolve three peaks [Hj]: two in the infrared region 
around 2000 cm" 1 (250 meV) and 5000 cm" 1 (0.6 eV), 
and one in the visible range around 1.6 x 10 4 cm -1 (2 
eV). Calculations in the framework of the t-J-Holstein 
model reproduce the middle peak at 0.6 eV IB"21| and, 
in more sophisticated approximations, also can distin- 
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guish the low energy peak at 2000 cm" 1 
ever, since this model projects out charge fluctuations, 
the highest peak at 2 eV is simply absent. Therefore, al- 
though the t-J-Holstein model is well suited to describe 
quasistationary states in the presence of a constant elec- 
tric field that is not too large [22j, since pumping by high- 
tion, and, of particular importance, the quantity probed energy pulses is associated with transitions into the up- 
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per Mott-Hubbard (or charge transfer) band, the HHM 
is the only model able to describe the time evolution of 
charge, magnetic, and lattice degrees of freedom in pump- 
probe experiments. We are aware of theoretical studies 
of pump-probe physics in models similar to the HHM in 
one dimension (ID) with classical [23| and quantum [24| 
treatment of phonons, and in 2D with a classical treat- 
ment of phonons [3, 0] ■ Here we present the first full 
quantum treatment of thermodynamic equilibrium and 
time-resolved optical spectroscopy of an electron-lattice 
system described by the HHM in 2D. Such calculations 
became feasible due to the recently developed double 
phonon cloud method (2(| which allows lattice degrees 
of freedom in large systems to be treated without uncon- 
trolled approximations. 

The HHM Hamiltonian coupled to an external classical 
vector potential is H = H t + Hy + Hepi, where 
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Hepi = w o ^ a\a% + 9^a ^(al + - n i)- (3) 
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Here t is the hopping amplitude, cj is the fermionic 
creation operator, is a unit vector along the axes of the 
lattice, a\ creates a phonon at site i with frequency ujo, 
and ni is the electron number operator. The EPI strength 
is defined by the dimcnsionless coupling constant A = 
g 2 uJo/4t. We choose model paremeters typical for high T c 
materials: t = 0.25 eV = 2000 cm" 1 and uj q = 0.2t = 400 
cm -1 . The value of the Hubbard repulsion U = Wt = 
2.5 eV yields low-energy physics very similar to that of 
the t-J model with J = At 2 /U = 800 cm -1 . The time- 
dependent potential vector is [Hj]: 



(■•■—■or 



A 



l e 



cos(w. 



pump \ 



to)), 



(4) 



and the OC, after the pulse, is given by [2j 
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<j(u,t) = —Z i e^ + ^ t (r\[j(t),j(0)}\r)dt 1 (5) 
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where |r) = Te'^o H (^) d ^ | r = 0), T is the time order- 
ing operator, j(t) is the current operator in the Heiscn- 
berg representation along one of the lattice direction 
axes, -q is a broadening factor taking into account ad- 
ditional dissipative processes, and M is the number of 
lattice sites. The |t) state is obtained through the Lanc- 
zos time propagation method [281 ]. In the following we 
use periodic boundary conditions on 2D lattices [29| with 
M = 10 for equilibrium and M = 8 for pump-probe 



calculations. The charge and spin degrees of freedom 
are evaluated exactly, whereas the quantum phonons are 
treated with the double phonon cloud method [26| . 

First we study the properties of the HHM in equilib- 
rium. The kinetic energy K and the average number 
of phonons N dressing the quasiparticle are plotted in 
Fig. [T] as a function of A, at half filling (HF) and at 
hole concentration 5 = 0.1 away from HF (one extra hole 
compared to HF in a system with M = 10 lattice sites). 
The strength of the EPI has a negligible effect at HF due 
to the uniform charge distribution, but in the doped sys- 
tem, the doped holes become lattice polarons at A c > 0.4 

We now turn to the OC of the HHM in equilibrium. 
Figure [21(a) compares the OC at zero (A = 0) and strong 
(A = 0.575) EPI in the half-filled system; Fig. \Sfb) 
presents the OC at zero and weak (A = 0.1) EPI at low 
hole concentration 5 = 0.1. The OC for strong EPI at 
half filling and at low hole concentration are compared in 
Fig.Efc). The OC of the HHM at half filling displays, in 
agreement with experiment (l5| . a gap with the peak of 
upper Hubbard band at u> ~ 8t m 16000 cm -1 . In con- 
trast to the t-J-Holstcin model, where this upper band is 
absent, the HHM can describe experiments using high- 
energy pump photons, ui pump > 8t. Similar to the ground 
state properties shown in Fig. [IJ the OC at HF weakly 
depends on A. The high-energy edge of the gap in OC 
is gradually filled by phonon-dressed excitations when A 
increases. 

On the contrary, the OC of the weakly doped system 
depends strongly on A as shown in Fig. [21(b). Upon chem- 
ical doping the gap is filled through a transfer of spectral 
weight towards lower energies. The most prominent fea- 
ture of the OC of the weakly doped system without EPI 
is a low energy peak at u ~ 2 J = 1600 cm -1 (solid green 
line in Fig. EJb)) which is well known from the pure t-J 
model. However, this narrow peak cannot explain the 
broad mid-infrared feature in the undcrdoped high T c 
compounds, because its energy is too low, the peak width 
is too narrow, and, for most model parameters, its dop- 
ing dependence is opposite that detected in experiments 
(see [li,|3il,|33l). Besides, the OC without EPI does not 
display the 2-peak structure of the experimental mid- 
infrared response [l9j . However, even a very small EPI 
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FIG. 1: (Color online) (a) Kinetic energy K and (b) average 
number of phonons JV versus A at HF and at 8 = 0.1. 
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(A = 0.1) induces a low-energy peak at u > ujo just above 
the phonon frequency (dashed black line in Fig.[2jb)). At 
the same time, the intensity of the 2 J peak is renormal- 
ized. Finally, the OC at larger values of EPI (A = 0.575) 
agrees well with the experimental positions of the OC 
peaks at energies 1500-2000 cm" 1 and 5000-6000 cm" 1 
(red solid curve in Fig. [2jc)) observed in strongly under- 
doped high T c compounds. The HHM model, like the t-J- 
Holstein model [19|, |34j , explains the doping dependence 



about 1/2 photon during the pulse with t p 
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of the spectra in the mid-infrared region. In addition, the 
HHM model describes the experimentally observed (lo| 
hardening of the large-energy «20000 cm" 1 peak associ- 
ated with the charge transfer gap with increasing doping 
(Fig. EJc)). Recently, it was argued that models with 
large Coulomb repulsion U, leading to Mott insulators, 
give OC with too little spectral weight at low energies, 
and a smaller value of U must be used in order to agree 
with experiments [Hj]. On the other hand, as shown in 
Fig. [2 EPI increases the low energy spectral weight of 
OC and, thus, large U Coulomb repulsion and EPI can 
provide another explanation of the effects discussed in 
Ref. [32| (see also j33|). Time-dependent pump-probe 
experiments can distinguish between these two scenar- 
ios. 

We now turn to the study of pump-probe dynamics, by 
solving the time dependent Schroedinger equation. For 
the chosen hopping t, the real time unit is H/t — 2.62 fs. 
We fix hujp Um p = 8.5i, above the peak of the charge trans- 
fer gap, and Aq = 0.2, corresponding to the absorption of 



and To = 5.24 fs. We study the dynamics in the first 250 
fs because our model is missing a heat reservoir where 
energy is dispensed in experiments during 500-1000 fs. 

The time evolution of the kinetic energy K, the av- 
erage value of Hubbard energy Ejj = (Hu), the lo- 
cal spin-spin correlation S = £^ „ Si ' Si+ M , the aver- 
age number of phonons N, and the average value of 
the contribution describing the charge-lattice coupling 

Ei = ( r lff w o J2i( a l + _ n -i)\ T ) are shown in Fig. [3] 
for three different values of EPI at half filling. Quan- 
tities associated with charge and spins (K, Ejj, and S) 
display quick changes during the first 10-20 fs (Fig. (3k., 
b, and c) followed by slower dynamics depending on the 
strength of the EPI. On the other hand, quantities asso- 
ciated with lattice degrees of freedom (N and £7) show 
dynamics on a time scale associated with the period of 
phonon oscillations T p h = 2tt/u!q ps 80 fs (Fig. [3ji) - Fur- 
thermore, after about two phonon oscillations, both at 
A = 0.1 and A = 0.5 a relaxed state is reached within a 
time around 150 — 200 fs. It is interesting to note that 
the spin correlation function S (Fig. [3b) at A / also 
displays oscillations at timescales similar to the phonon 
period T p h 80 fs, highlighting the strong coupling of 
spin and lattice degrees of freedom. Notably, after the 
laser pulse, slower charge dynamics are observed only in 
the presence of a significant EPI (see inset Fig. |3]b) : at 
A = 0, Eu diplays oscillations [IH around an essentially 
constant value, whereas at A = 0.5, a relaxation is ob- 
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FIG. 2: (Color online) OC at equilibrium (in units of 00 
e 2 /ha) for different EPI values, at HF and at 8 = 0.1. 
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FIG. 3: (Color online) Time evolution of: (a) Kinetic energy 
K; (b) and inset in (b): U interaction contribution Eu = 
(Hu); (c) and inset in (c): local spin-spin correlation; (d): 
average number of phonons iV; inset in (d): EPI energy at 
A = (dotted green line), A = 0.1 (solid black line), and 
A = 0.5 (dashed red line). 
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served. The relevance of the EPI is also indicated by 
strong A dependence of K after the pumping: the kinetic 
energy gain, that is A independent in the ground state, 
quickly reduces by increasing A(35| . 

The oscillation behavior with characteristic phonon 
frequenciespredicted here, has not been observed in ex- 



pcrimcnt [5= 



To our mind, the difference lies 



in the pump-width to phonon period ratio T pum p/T p h 
which is small T pump /T p h <C 1 in our calculation and 
large T pump /T p h > 1 in experiments 043 l^-HH- This 
interpretation is supported by numerous pump-probe ex- 
periments in charge-transfer organic compounds fl2- 14 1. 
where the relation T pump /T p h <C 1 holds, and the os- 
cillations of physical quantities with phonon period are 
observed. Moreover, experiments on high T c materials 
with record ultra-short T pump = 12 fs pulse [1] do ob- 
serve oscillations with the frequencies of bond-bending 
and apical oxygen phonons. In addition, the wavelet 
transformation of time dynamics of the OC shows that 
the energy transfer into the lattice system is associated 
with a time around r tr ~ 150 fs, which is similar to that 
of phonon oscillations Hence, interplay of similar 
periods r tr ~ T p h may result in beat behavior of lat- 
tice properties as observed in our calculations at A = 0.5 

(Fig. EH). 

Fig|4] shows the OC after the pulse at HF with (a) 
zero, (b) weak, and (c) strong EPI. The gap present in 
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FIG. 4: (Color online) Non-equilibrium OC at three different 
times after the pump pulse at (a) A = 0, (b) A = 0.1 and (c) 
A = 0.5. Inset shows zoom of bimagnon feature in the OC at 
HF and at equilibrium. 



the undoped system is filled after the pulse, but in a way 
that is strongly dependent on the EPI strength. The 
non-equilibrium OC without EPI shows almost time in- 
dependent Drude behavior (Fig.@Ji)- On the other hand, 
for weak EPI (A = 0.1) as shown in Fig. gb, the OC 
oscillates and contains two main peaks at the same en- 
ergies as observed in the weakly doped system at weak 
EPI (Fig. Eb)- These peaks correspond to 2 J and ujq, 
which are related to charge scatterings with magnons 
and phonons, respectively. These two features, associ- 
ated with spin and lattice degrees of freedom, become 
visible at very short times after the pump pulse. This 
result agrees with results obtained within the Holstein 
model [38|. Finally, the non-equilibrium OC at A = 0.5 
has no features in common with the OC at equilibrium 
(compare Fig.0h and Fig. [5];): there is no dominant peak 
at 5000 cm" 1 or at 2000 cm" 1 . Instead, a low energy 
peak with phonon sidebands is observed indicating that 
the lattice polaron is not yet formed. Indeed polaron 
formation at intermediate EPI requires long times [38|. 
due to the high potential barrier between free and self- 
trapped states. 

At zero (weak) EPI, a feature is clearly (tentatively) 
revealed in the non-equilibrium OC at u ~ 3 J, due to the 
transition into the bimagnon state. This excitation has 
been observed by Perkins et al. j3(| in undoped copper 
oxides and, theoretically, can be seen in the equilibrium 
OC only in the presence of EPI [37| (see inset in Fig. 
for equilibrium OC at different EPIs). However, out of 
equilibrium, due to lower symmetry, such a transition is 
allowed even at A = 0. 

In conclusion, we studied spectral properties of under- 
doped high T c materials (2D Hubbard-Holstein model 
with quantum phonons) in equilibrium and during the 
first 250 fs following excitation by a laser pulse. In partic- 
ular, we found phonon-govcrned oscillations that can be 
observed in pump-probe experiments provided that the 
pulse time T pump is shorter or at least comparable with 
the typical phonon times T p h = 2it/uj . We showed that 
the analysis of OC in and out of equilibrium in the mid- 
infrared range can reveal the EPI strength in cuprates. 
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